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ABSTRACT 

Thermal barrier coating life models developed 
under the NASA Lewis Research Center's Hot Section 
Technology (HOST) program are summarized. An initial 
laboratory model and three design-capable models are 
discussed. Current understanding of coating failure 
mechanisms are also summarized. 

INTRODUCTION 

Thermal barrier coatings are being developed for 
protecting air-cooled turbine blades and vanes In air- 
craft gas turbine engines. The current state-of-the-art 
coating system consists of about 0.25 mm of a zlrconla- 
yttrla ceramic over 0.13 mm of an MCrAlY alloy bond 
coat. Both layers are applied by plasma spraying. The 
benefits arise from the Insulation provided by the 
ceramic layer. This Insulation allows higher gas tem- 
peratures, lower component temperatures, reduced cool- 
ing air requirements, moderation of thermal transients, 
and/or a decrease In the severity of hot spots. This 
yields Improvements In performance, efficiency, and com- 
ponent durability. Future engine designs are expected 
to rely heavily on thermal barrier coatings. Thus life 
models are required to assess the risks associated with 
any given design and to insure that these coatings can 
be exploited fully. Further details may be found In 
Miller (1987), DeMasl et al., (1988), Strangman et al., 
(1987), and Hlllery (1987). 

The NASA thermal barrier coating life model devel- 
opment program consisted Initially of an In-house pro- 
gram designed to improve understanding and to develop a 
model suitable for treating laboratory life data 
(Miller, 1987). This work was then extended via three 
contracts under the HOST program to the development of 
design-capable models (DeMasl et al., 1988; Strangman 
et al., 1987; Hlllery et al., 1987). These contracts 
were devised determine thermomechanical properties, to 
analyze coating stresses and strains, and to develop 
life models. Phase I of each contract has now been suc- 
cessfully completed, and the results will be summarized 
In this paper. 


COATING FAILURE MECHANISMS 

A basic understanding of coating failure mechanisms 
Is a prerequisite to the development of life prediction 
models. Failure mechanisms In gas turbine engines and 
In laboratory simulations have been discussed In detail 
elsewhere (e.g.. Miller, 1987; DeMasl et al., 1988; 
Strangman et al., 1987, Hlllery et al., 1987). There 
Is now general agreement that these coatings fail pri- 
marily as a result of stresses Induced by the thermal 
expansion mismatch between the ceramic and metallic 
layers, and that these stresses are greatly influenced 
by tlme-at-temperature processes such as oxidation and 
possibly sintering. The stress state In the ceramic 
layer which leads to crack propagation and eventual 
spalling Is one of biaxial compression In the plane of 
the coating and radial tension. These stresses are fur- 
ther complicated by the wavy and Irregular interface 
between the ceramic and metallic layers. In fact, 
HOST-sponsored calculations indicate that the radial 
stresses above a wavy interface may actually alternate 
between regions of compression and tension as illus- 
trated in Fig. 1 (Chang et al., 1987). Figure 2 
Illustrates that the behavior of plasma sprayed 
zlrconla-based thermal barrier coatings differs signifi- 
cantly from the behavior of conventional ceramics. This 
behavior, which Is believed to result from the splat 
structure, includes very low thermal conductivity and 
very high compressive strain tolerance. In-plane ten- 
sile strain tolerance of the coating system is also very 
high because such loading may lead to segmentation 
cracking in the ceramic with no degradation to the 
attachment strength. Plasma sprayed zirconla-yttrla 
also exhibits creep-like behavior, presumably as a 
result of sliding at the splat boundaries, and fatigue- 
like behavior, presumably as a result of slow crack 
growth. Experimental evidence of slow crack growth (or 
microcrack link up), creep, and fatigue are presented 
In DeMasl et al., (1988). 

INITIAL LABORATORY MODEL DEVELOPMENT 

A preliminary life prediction model has been 
described (Miller, 1987; Miller, 1984; Miller et al . , 
1984). This model assumed that the complex state of 
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stress and strain imposed on the coating system by the 
thermal loads could be expressed in terms of a single 
parameter. This parameter was labelled e r — which was 
taken to be the radial component of the thermal expan- 
sion mismatch strain. Next it was assumed that the 
time-at-temperature effects could be treated in terms 
of oxidation alone and that oxidation could be charac- 
terized by the weight gain at the conclusion of eacn 
cycle wn. Then, weight gain and strain were related 
using either of two alternate approaches. In the first 
case, depicted in Fig. 3(a), an oxidized coating is 
assumed to behave as if an effective strain e e is 
increasing. At zero weight gain this effective strain 
equals the radial strain e r . At a critical weight gain 
w c — defined as the weight gain required to fail the 
coating in a single cycle — the effective strain equals 
a failure strain ef. This leads to the expression 



where the exponent m has been added to allow the curve 
in Fig. 3(a) to be nonlinear. The alternate assumption 
(Miller, 1987) is to assume that the failure strain 
degrades from an initial value efo to a final value 
equal to e r . This case, illustrated in Fig. 3(b) leads 
to the expression 



.xhere Nfjyj is the apDarent number of cycles remaining 
after cycle N and weight gain 

Figure 4 illustrates the fits obtained using 
expression 4 and applying it to life aata collected at 
i 1 00 °C for three different cycle lengths. It should 
be mentioned that the set of parameters given in the 
figure are not unique. Numerous other sets provide 
equally good fits. For example raising the assumed 
value of b wnile lowering the strain ratio produces 
an equally ;ood fit. Also, the life data can be fit 
equally we i 1 using expressions 4 or 5. 

DESIGN-CAPABLE LIFE MODELING 


While the above model represented a first step it 
was net in a form which would be of use to an engine 
designer, "herefore three contracts were instituted 
under the HOST program which were aimed at the develop- 
ment of design-capable models. 

Pratt & Whitney Aircraft (DeMasi et al . , 1988), 
along with subcontractor Southwest Research Institute, 
developed a fatigue-based coating life model which uses 
Miner’s Law (expression 7) along with expression 6 
rewri tten as 



wnere Ae} 


is the inelastic strain range defined by 
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(9) 


Cracks in the ceramic layer may be assumed to grow 
according to a crack growth law of the form 


dA 

dN 


AeV 
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(3) 


where da/dN is the incremental crack growth per cycle, 
A is a constant, b and c are exponents related to 
the subcritical crack growth exponent, and a is the 
crack length. The model resulting from expression 1 is 



N.= l 

and the alternative model resulting from expression 2 is 



N = 1 


These models may also be derived from the familiar 
fatigue expression (Miller et al., 1984; Manson, 1966). 



and Miner's Law 



N = 1 


( 6 ) 


(7) 


The term A(a AT) in the above expression is the ther- 
mal expansion mismatch strain (which was expressed in 
terms of e r in the previous section), Ach is the 
strain resulting from the heating transient, Ae c is 
the strain resulting from the cooling transient, and 
a vs /E is the elastic strain at yielding. The assumed 
relationship between oxidation and strain, analogous to 
expression 2, was 



where oxidation has been expressed in terms of the 
oxide layer thickness $ rather than the specific 
weight gain w. The inelastic strain range was calcu- 
lated using finite element techniques which employed a 
time dependent inelastic model developed by Walker 
(1983). Figure 5 shows an example of the use of this 
model to calculate compressive and tensile strains which 
may be compared with experimental data. In Fig. 6 the 
ceramic stress-strain behavior is calculated for a sin- 
gle cycle. This figure displays the large amount of 
reversed inelastic strain produced by thermal cycling. 
Figure 7 shows a plot of observed versus calculated 
lives for a wide range of test conditions. As shown in 
the figure the model is accurate to plus or minus a fac- 
tor of 3, which is considered adequate. 

The model developed by the Garrett Turbine Engine 
Company (Strangman et al . f 1987) may be expressed as 
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Equation 11 is expressed schematically in Fig. 8 which 
shows that the model considers bond coat oxidation, zir- 
conia toughness reduction, and damage due to molten salt 
deposits. The model is driven by the thermal analysis 
of the component of interest for its anticipated mis- 
sion. The left side of the denominator in expression 11 
as determined from test data calibrations is 

0.25 

(t + 0.181) MTBREF 

(exp[-0.015(T ♦ 273) + C^}' 1 + (exp[-0.041( -*-273) + C^}” 1 

( 12 ) 

where MTBREF is a mul ti temperature burner rig experi- 
ence factor which forces predictions and experiment into 
agreement. The right side of the denominator in expres- 
sion 11 is calculated using a Garrett-developed model 
(Strangman, 1984; Strangman et al. f 1987). In practice, 
the model is driven by thermal analysis of the component 
of interest. An example of the application of the ther- 
mal barrier coating life model to laboratory test data 
is shown in Fig. 9, and mission analysis predictions are 
shown in Fig. 10. 

The approach used by the General Electric Company 
(Hill ery et al . , 1987) employed time-dependent, nonlin- 
ear finite element modeling of the stresses and strains 
present in the thermal barrier coating system, followed 
by the correlation of these stresses and strains with 
test lives. The life model developed using this 
approach may be expressed as 

Ac rz + 0.4 Ae R = 0.121 N f ' 0-486 (13) 

where Agrz is the shear strain range, Acr is the nor- 
mal strain range, and Nf is the number of cycles to 
failure. The above model is the only one to consider 
failure induced by edges and hence is the only one to 
consider shear strain. Expression 13 is Illustrated 
graphically in Fig. 11 

CONCLUDING REMARKS 

In conclusion, the materials and structural aspects 
of thermal barrier coatings have been successfully inte- 
grated under the NASA HOST program to produce models 
which may now or in the near future be used in design. 
Efforts on this program continue at Pratt & Whitney Air- 
craft where their model is being extended to the life 
prediction of physical vapor deposited thermal barrier 
coatings . 

While the HOST program has been quite successful it 
should also be noted that many new and unanswered ques- 
tions have been raised by this work. For example, the 
effects of creep and inelasticity in both the ceramic 
and bond coat layers are poorly understood. The role of 


shearing stresses, including the role that shearing at 
an edge may play in reducing the fatigue exponent, is 
not well understood. The detailed mechanism by which 
oxidation controls coating system life is not well 
understood either. Also, it is not known whether the 
assumption of a smooth interface, commonly employed to 
simplify finite element analyses can lead to inaccurate 
or even misleading results. Other areas of uncertainty 
involve the importance of sintering at high temperatures 
and hot corrosion at relatively low temperatures. 
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FIGURE 1. - SCHEMATIC REPRESENTATION OF CALCULATED RADIAL 
THERMAL EXPANSION MISMATCH STRESS ABOVE A WAVY INTER- 
FACE. 
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FIGURE 2. - SCHEMATIC REPRESENTATION OF THERMAL MECHANICAL 
PROPERTIES RESULTING FROM COATING SPLAT STRUCTURE. 
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FAILURE STRAIN- £ f EFFECTIVE STRAIN, 



(A) ASSUMED RELATIONSHIP BETWEEN EFFECTIVE STRAIN 
AND OXIDATIVE HEIGHT GAIN. 



(B) ASSUMED ALTERNATE RELATIONSHIP BETWtEN FAILURE 
STRAIN AND OXIDATIVE HEIGHT GAIN. 

FIGURE 3. - ASSUMED HEIGHT GAIN/STRAIN RELATIONSHIPS. 
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FIGURE A. - COMPARISON OF CALCULATED AND MODELED LIFE AS 
A FUNCTION OF HEATING CYCLE DURATION ACCORDING TO NASA 
LABORATORY MODEL. 



FIGURE 5. - CERAMIC BEHAVIOR MODELED WITH WALKER EQUATION 
ACCORDING TO THE PRATT t WHITNEY MODEL. 
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FIGURE 6. - SCHEMATIC OF STRAINS CALCULATED FOR A TYPICAL 
CYCLE ACCORDING TO THE PRATT & WHITNEY MODEL. 



FIGURE 8. - SCHEMATIC OF PROCESSES CONSIDERED IN THE 
GARRETT MODEL. 



OBSERVED CYCLE LIFE 

FIGURE 7. - COMPARISON OF CALCULATED AND 
MODELED LIVES ACCORDING TO THE PRATT & 
WHITNEY MODEL. 
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FIGURE 9. - CONSERVATIVE FIT OF LABORATORY LIFE DATA 
TO THE GARRETT MODEL. 
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FIGURE 10. - MISSION ANALYSIS PREDICTIONS 
BY THE GARRETT MODEL. 
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FIGURE 11. - CORRELATION BETWEEN CALCULATED STRAIN RELA- 
TIONSHIP AND EXPERIMENTAL CYCLES TO FAILURE ACCORDING 
TO THE GENERAL ELECTRIC MODEL. 
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